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SECTION  I 


INTRODUCTION 

High  performance  gas  turbines  require  high  turbine  inlet  temperatures 
while  the  skin  of  exposed  components  has  to  be  maintained  at  a  lower  tem¬ 
perature  for  long  life  expectancy. 

This  lower  skin  temperature  can  be  obtained  by  a  combined  use  of 
internal  convection  cooling  and  of  an  external  thermal  barrier  in  the  form 
of  a  film  of  cool  air,  covering  the  skin  of  exposed  components.  The  design 
of  such  turbines  requires  detailed  information  on  the  heat  transfer  distri¬ 
bution  and  aerodynamics  of  highly-loaded  stages,  and  on  the  effectiveness 
of  the  film  thermal  barrier.  This  film  can  be  obtained  by  blowing  air 
through  a  porous  surface,  a  slot,  or  an  array  of  holes.  The  first  tech¬ 
nique,  transpiration  cooling,  requires  relatively  small  cooling  flows.  The 
momentum  of  the  blown  air  is  negligible  compared  to  the  main  stream  momen¬ 
tum,  and  the  boundary  layer  theory  can  be  applied  to  determine  the  thermal 
effectiveness  of  the  film.  However,  the  technological  implementation  of 
transpiration  cooling  has  met  severe  difficulties,  and  the  second  technique 
of  film  cooling  has  found  many  applications  in  advanced  gas  turbines. 

Since  the  air  is  injected  through  a  limited  number  of  orifices,  often  at  an 
angle  to  the  wall,  momentum  of  the  injected  flow  is  not  negligible  compared 
to  the  main  flow,  and  the  boundary  layer  theory  cannot  be  applied  to  this 
case.  New  analytical  and  experimental  information  is  therefore  needed  to 
determine  the  effectiveness  of  the  film  thermal  barrier  obtained  by  a  wide 
variety  of  injection  configurations. 
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The  critical  information  needed  for  film  cooling  deals  with  coolant 
penetration  into  the  main  flow  and  with  the  generation  of  streamwise  vor- 
ticity  in  the  vicinity  of  the  injection  point.  Further  downstream, 
information  is  needed  to  understand  the  coupled  turbulent  and  vortical 
mixing  of  the  coolant  with  the  mainstream. 

Engine  performance  optimization  depends  critically  on  the  selection  of 
an  effective  cooling  configuration  and  a  detailed  understanding  of  heat 
transfer  to  turbine  components.  A  systematic  study  of  relevant  heat 
transfer  and  film  cooling  effectiveness  is  described  in  this  report,  which 
covers:  O)  analytical  and  experimental  studies  of  a  wide  variety  of  film 
cooling  configurations,  (2)  experimental  and  analytical  predictions  of  heat 
transfer  over  highly-loaded  profiles,  and  (3)  film  cooling  effectiveness 
over  a  subsonic  profile.  Short  duration  studies  using  a  shock  tunnel  and 
hot  blow-down  facilities^  provide  a  practical  and  flexible  means  to 
acquire  the  necessary  experimental  information  while  providing  good 
modeling  of  the  flow  and  temperature  fields.  These  short  duration  studies 
are  validated  by  long  duration  high  temperature  tests  of  an  air  cooled  flat 
plate. 
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2.1  Adiabatic  and  Isothermal  Wall  Heat  Transfer 

Heat  transfer  measurements  are  conventionally  referred  to  as  adiabatic 
conditions.  The  determination  of  heat  transfer  rate  requires  the  knowledge 
of  film  coefficient  and  of  an  adiabatic  wall  temperature  (Tgc)).  The  heat 
transfer  coefficient  is  determined  for  the  case  of  no  injection  and  a  well- 
developed  boundary  layer.  An  adiabatic  wall  effectiveness  is  defined  as 


T  .  -  T 
ad  os 

nad  5  T  -  T 

C  09 

where  T^  and  TQ  are  the  main  stream  and  coolant  temperatures.  This 
adiabatic  approach  is  correct  when  the  ratio  of  the  main  stream  to  film 

temperatures  is  close  to  unity.'  In  advanced  gas  turbines,  heat  diffusion 
brings  metal  temperatures  closer  to  an  isothermal  temperature  distribution 
rather  than  an  adiabatic  one.  In  the  presence  of  angular  injection  and 
high  main  flow  to  coolant  temperature  ratios,  and  in  close  vicinity  to  the 
injection  point  where  large  gradients  of  flow  properties  exist,  the  vali¬ 
dity  of  using  adiabatic  wall  information  to  determine  heat  loads  to  the 
cooled  surfaces  is  in  doubt,  thus  requiring  the  investigation  of  heat 
transfer  phenomena  by  directly  measuring  heat  transfer  under  correctly 
modeled  turbine  conditions.  This  need  provides  the  justification  to  refer 
the  data  to  isothermal  conditions.  Under  these  conditions,  the  heat  trans¬ 
fer  rate  is 


Q/A  =  h(T  -  T  ) 

W  00 


where  Tw  is  the  temperature  of  the  isothermal  wall  and  the  isothermal 
effectiveness  is  defined 
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n,„  =  1  -  Q./Q. 
iso  c  o 

where  Qc  is  the  heat  transfer  with  cooling  and  Qq  is  the  heat  transfer  at 
the  same  condition  without  cooling.  The  isothermal  effectiveness 
represents  the  capacity  of  the  thermal  barrier  provided  by  the  film  to 
reduce  the  heat  transfer. 

Adiabatic  and  isothermal  effectivenesses  are  related  by 


Tc  -  T 

”130  =  1  -  Vho  <’  *  "ad  T-TT> 

00  W 

when  h„  and  h  are  the  film  coefficients  with  and  without  cooling, 
co 

2.2  Experimental  Techniques 

Unsteady  heat  transfer  measurements  under  steady  flow  conditions  pro¬ 
vide  the  advantages  of  relatively  simple  heat  transfer  measurements  of 
short  duration,  and  of  flexibility  in  modeling  actual  gas  turbine  operating 
conditions.  To  model  gas  turbines,  one  needs  a  high  enthalpy  flow  for  the 
mainstream  and  a  low  enthalpy  flow  for  the  coolant.  By  scaling  the  temper¬ 
atures  of  these  flows  and  the  skin  temperature  from  the  actual  turbine  con¬ 
ditions  and  matching  Reynolds,  Mach  and  Prandtl  numbers,  a  fairly  represen¬ 
tative  model  can  be  obtained  for  the  heat  transfer  and  cooling  problems. 

The  exDerimental  technique  used  is  to  establish  the  high  and  low  tem¬ 
perature  flows  concurrently  in  a  stepwise  fashion  for  a  duration  short 
enough  to  render  wall  temperature  changes  minimal  so  that  heat  fluxes  can 


be  measured  on  a  quasi-isothermal  surface. 

The  short  duration  necessitates  the  use  of  fast  response  instrumenta¬ 
tion  such  as  thin-film  gauges,  calorimetric  gauges,  and  Kistler  pressure 
transducers. 
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A  shock  tunnel  and  a  blow-down  facility^  are  used  for  these  tests. 
The  shock  tunnel  is  used  to  obtain  basic  data  on  film  cooling  effectiveness 
over  flat  and  curved  plates.  The  blow-down  facility  can  be  used  either  in 
conjunction  with  a  cascade  or  with  a  rotating  turbine  rig  to  determine 
distribution  of  heat  transfer  rates  and  film  cooling  effectiveness  over 
blade  profiles,  end  walls,  and  turbine  shrouds. 

2.3.1  The  shock  tunnel^  (Figure  1)  is  attached  to  a  3"  shock  tube 

through  an  accelerating  nozzle.  At  the  downstream  end  of  the  plexiglass, 

constant  area  (2  in.  x  7/8  in.)  tunnel,  a  choked  nozzle  controls  the  Mach 

number.  The  steady  flow  test  time  is  approximately  10  ms,  whereas  the  time 

to  establish  flow  is  1  ms.  For  those  short  durations,  tests  of  the  order 

of  10  ms,  the  plexiglass  walls  have  large  enough  thermal  inertia  to  limit 

their  temperature  rise  to  less  than  7°F,  thus  not  offsetting  the  isothermal 

wall  boundary.  Various  configurations  of  coolant  holes  (or  slots)  on  the 

tunnel  walls  are  connected  to  a  coolant  flow  supply  plenum.  The  location 

5 

of  the  holes  (or  slots)  provides  a  local  flow  Reynolds  of  3  x  10  .  The 
coolant  flow  is  prepurified  nitrogen,  which  is  cooled  by  passing  it  through 
a  copper  coil  immersed  in  a  bath  of  liquid  nitrogen.  The  introduction  of 
the  coolant  is  synchronized  with  the  rupture  of  the  shock  tube  diaphragm. 
Arrays  of  thin-film  gauges  are  located  downstream  of  the  injection  point  to 
determine  the  heat  transfer  rate  as  a  function  of  distance. 


2.3.2  The  Hot  Blow-down  Facility^  (Figure  2)  consists  of  air  supply 
system  storing  400  lbs  of  air  at  2400  psi ,  dome-loading  type  pressure  reg¬ 
ulators,  a  6  ft  long  and  1.5  ft  diameter  pebble  bed  heater  filled  with  3500 
lbs  of  stainless  steel  balls  and  a  filter  trap.  The  facility  has  two 
branches,  each  consisting  of  a  diaphragm  section,  a  metering  throat,  a 
pressure  reducing  section  utilizing  a  supersonic  nozzle  and  diffuser,  and  a 
transitional  piece  leading  in  one  branch  to  the  turbine  rig  and  in  the 
other  branch  to  the  cascade.  The  tests  are  typically  less  than  1  s, 
whereas  steady  flow  is  established  after  0.15s  -  0.25s.  The  heat  transfer 
rate  is  measured  by  calorimetric  gauges,  a  thin  aluminum  or  copper  disc 
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where  n  is  a  number  like  0.5  to  0.7,  one  obtains 


P  L  T 

op  .  op  =  op 

PR  LR  TR 
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=  (S) 


n  +  1/2 


The  scaling  implies  that  if  room  temperature  is  used  to  model  turbine 

components  temperatures  which  are  as  high  as  1650°F,  the  temperature 

scaling  factor  would  be  4.  If  the  length  scale  is  kept  the  same  and 

1*14 

n  =  0.64,  the  pressure  would  be  reduced  by  the  ratio,  Pop/Pr  =  s  =  **.36. 
For  the  tests  with  film  cooling,  the  temperature  of  the  coolant  is  reduced  py 


the  ratio  S  =  4. 


SECTION  3 

3 .  Film  Cooling  Effectiveness 

As  mentioned  before,  engine  designs  geared  to  long  life  do  not  pre¬ 
sently  use  transpiration  cooling,  seldom  use  injection  parallel  to  the 
mainstream  and,  in  most  cases,  use  injection  at  an  angle  to  the  mainstream. 
As  a  reference,  however,  the  state  of  understanding  of  parallel  injection 
will  be  reviewed  briefly. 


3. 1  Parallel  Injection 

When  the  film  is  injected  parallel  to  the  main  flow,  turbulent  mixing 

dominates  the  interaction  between  the  two  streams.  Downstream,  the  mixed 

layer  is  dominated  by  the  wall  shear  stress,  therefore,  by  suitable  balance 

(2) 

of  both  energy  and  mass  for  the  mixed  streams,  Stoliery  and  El  Ehwany  , 

(3) 

and  Leontev  and  Kutateladze  have  derived  the  following  relation: 
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applicable  for  x>>s,  where  s  is  the  slot  height,  with  A  =  3. 03  in  reference 
2,  and  A  =  3.1  in  reference  3,  where  Rgc  is  the  Reynolds  number  based  on 
coolant  conditions. 

This  theoretical  relationship  compares  favorably  with  the  test 
results^^^  for  data  obtained  far  downstream,  but  the  scatter  becomes 
appreciable  in  the  region  closer  to  the  injection  point.  In  the  vicinity 
of  the  injection  point,  the  mixing  is  dominated  by  the  local  vorticity 
field  which  is  not  included  in  the  above  theoretical  formulation. 
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3.2  Film  Cooling  with  Angular  Injection 

When  a  film  is  injected  at  an  angle  to  the  mainstream  from  a  well  to 
provide  a  thermal  barrier  over  a  section  of  wall,  it  undergoes  two  deter¬ 
mining  processes:  a  viscous  interaction  and  a  pressure  interaction.  The 
viscous  interaction  results  in  the  turbulent  mixing  of  the  two  streams. 
The  pressure  .nteraction  exists  since  the  momentum  of  the  injected  flow 
normal  to  the  wall  is  commensurable  with  the  main  stream  momentum.  The 
pressure  gradient  normal  to  the  wall  near  the  injection  point  is  signi¬ 
ficant  and  boundary  layer  theories  are  no  longer  applicable. 

The  pressure  interaction  determines  the  penetration  of  the  coolant 
into  the  boundary  layer  and  in  the  mainstream,  and  also  the  formation  of 
streamwise  vorticity.  The  viscous  interaction  determines  the  rate  of 
entrainment  of  the  mainstream  fluid  by  the  coolant  and  the  rate  of  mixing 
of  the  two  streams. 

Further,  the  pressure  interaction  for  injection  parameters  commonly 
used  leads  to  a  separation  of  the  coolant  and  boundary  layer  at  the 
injection  point  and  leads  to  the  formation  of  a  separation  bubble  just 
downstream  of  the  injection  point.  The  evidence  for  the  existence  of  the 
boundary  layer  separation  at  the  injection  point  and  the  formation  of  a 
separation  bubble  are  several: 


3.2.1  This  separation  is  theoretically  expected  when  the  criterion 

(7) 

determined  by  Kutateladze  and  Leontev  is  applied  to  the  conditions  of 
film  cooling  with  angular  injection.  This  criterion  indicates  that 
separation  will  occur  when  the  normal  mass  injection  parameter 
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where  Cfo  is  the  skin  friction  coefficient  in  the  absence  of  coolant,  and  y 
refers  to  the  direction  normal  to  the  wall.  For  a  flow  Reynolds  number  at 
the  injection  point  of  Rg  =  3  x  10'’,  C^q  =  2.5  x  10  and  this  criterion 
would  indicate  that  separation  is  likely  to  occur  when  my  »  5  x  10“3, 
i.e.,  for  a  typical  angle  of  injection  of  30°,  separation  would  occur  for  a 
value  of  the  total  mass  injection  parameter. 


o„  U 

c  c 

m  =  >  0.1 

p*p*  ~ 

03  CO 


Practical  values  of  the  injection  parameter  are  generally  appreciably 
larger  than  m  =  0.1. 


3.2.2  The  existence  of  a  separation  bubble  downstream  of  30  slot  was 

demonstrated  by  dye  and  hydrogen  bubble  visualization  techniques  in  a  water 
(  ft  ) 

tunnel  .  The  study  indicated  that  the  length  of  the  separation  bubble 
was  proportional  to  the  momentum  ratio 
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The  flow  visualization  also  indicated  the  appearaice  of  streamwise  vortices 
having  consecutive  opposite  directions  of  rotation.  This  vorticity  results 
from  the  coupling  between  the  injected  flow  and  its  streamwise  velocity. 

The  vortical  motion  has  the  effect  of  enhancing  the  mixing  of  the  coolant 
with  the  mainstream  as  evidenced  by  the  presence  of  fluid  elements  in  the 
separation  bubble.  The  study  also  indicated  that  the  vortices  persisted 
far  downstream. 
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3.2.3  Measurements  of  heat  transfer  coefficient  (Figure  3)  downstream 
of  a  30°  injection  slot  on  a  curved  wall  indicated  that  the  heat  transfer 
coefficient  goes  through  a  maximum  close  to  the  injection  for  a  given  value 
of  the  injection  parameter.  This  maximum  can  be  explained  by  the  existence 
of  the  separation  bubble  in  which  hot  gas  is  being  entrained. 

3.2.4  Measurements  of  heat  transfer  rate  at  the  trailing  edge  of  a 

(q) 

plate  indicate  that  the  isothermal  effectiveness  measured  at  the 

trailing  edge  can  be  severely  reduced  (Figure  4)  by  locating  the  coolant 
injection  slot  at  a  distance  shorter  than  the  reattachment  point  of  the 
separation  bubble  (Figure  5).  The  separation  bubble  is  observed  in  the 
range  of  mass  injection  up  to  values  of  m,  for  which  complete  separation  of 
the  boundary  layer  occurs. 

In  the  case  of  injection  through  inclined  holes,  the  separation  down¬ 
stream  of  the  injection  point  is  mitigated  by  the  formation  of  counter 
rotation  vortices  which  persist  far  downstream  and  contribute  to  the 
downstream  mixing  oi  the  two  streams. 

As  it  will  be  described  below,  streamwise  vorticity  is  also  generated 
when  slot  injection  is  used. 

3-3  Film-Cooling  Effectiveness  with  Angular  Injection 

A  comprehensive  experimental  program^®^ 1 1  ^ ^  has  been  undertaken 
to  investigate  different  film  cooling  configurations  of  slots,  single  and 
double  line  holes  for  different  streamwise  and  crosswise  angular  injection. 

Figure  6a  is  a  representative  layout  of  the  flat  plate  used  in  the 
shock  tunnel  for  the  evaluation  of  film-cooling  effectiveness.  The  injec¬ 
tion  scheme  uses  a  double  line  of  holes  which  make  an  angle  c^,  with  the 
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streamwise  direction  and  a ^  with  the  crosswise  direction.  Thin-film  heat 
transfer  gauges  are  located  at  different  distances  corresponding  to  dif¬ 
ferent  values  of  the  ratio  of  the  streamwise  distance  (x)  to  the  hole  dia¬ 
meter.  These  gauges  provide  an  average  local  heat  transfer  rate,  whereas 
the  gauges  located  at  x/D  =  27  are  used  to  measure  the  crosswise  variation 
of  effectiveness.  Most  of  the  tests  were  conducted  at  a  stagnation  temper- 
atue  of  1000°R  and  10  psia  stagnation  pressure,  which  simulated  the  inlet 
conditions  ilOOO°R  and  175  psia  of  an  advanced  gas  turbine.  The  Mach  number 
was  typically  M  =  0.5  and  the  coolant  temperatures  were  220°R,  310°R,  400°R 
and  500°R.  The  mass  injection  parameter  was  varied  between  0.1  to  1.6. 

The  streamwise  angles  used  were  10°,  20°  and  30c  and  the  crosswise  angles 
used  were  0°,  30°,  50°,  70°,  and  90°. 

The  influence  of  the  different  configuration  parameters  can  be 
assessed . 

3.3.1  Effect  of  the  Crosswise  Angle  for  Holes 

The  cross-stream  injection  angle  does  not  effect  the  film  cooling  ef¬ 
fectiveness  over  a  flat  plate  to  any  great  degree.  Its  influence  is 
limited  to  the  area  just  downstream  of  the  injection  point  as  indicated  by 
the  film  effectiveness  footprints'11^,  e.g.,  (Figure  6b),  measured  for 
-  30°  and  for  different  values  of  the  mass  injection  parameter.  In  the 
region  close  to  the  hole,  the  coverage  increases  with  the  crosswise  width 
of  the  hole. 

3.3.2  Influence  of  the  Distribution  of  Holes 

Data  of  film  effectiveness  were  correlated  by  using  the  concept  of 
equivalent  slot  width,  which  is  the  width  of  a  slot  which  would  pass  the 
same  amount  injection  mass  flow  as  the  single  or  double  line  of  holes. 
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A*  =  K  _1,35  R  '*25  x/s 
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where  K  is  the  square  root  of  the  momentum  ratio  and  R  is  the  Reynolds 
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number  based  on  the  injection  conditions. 


3.3.3  Slot  versus  Hole  Injection 

As  seen  in  Figure  8,  the  film  cooling  effectiveness  of  hole  injection 
is  sharply  lower  than  injection  through  a  slot,  corrected  for  equivalent 
coolant  mass  flow  introduced  to  form  the  film.  By  observing  that  the 

mixing  area  for  a  circular  jet  is  w  times  larger  per  unit  length  of  3lot, 

» 

the  correlation  parameter  AQ  was  corrected  by  this  amount  for  data  related 
to  holes.  In  other  words  irX  =  1.0  for  a  slot  and  *  =  u  for  hole  injec¬ 
tion.  With  this  change  in  the  correlation  parameter,  the  hole  data  aligned 
with  the  slot  data  quite  well  (Figure  9).  The  best  overall  agreement  is, 

x 

however,  obtained  when  the  factor  f  is  raised  to  the  power  1.41. 


3.3.4  Overall  Correlation 

With  the  incorporation  of  this  parameter  and  the  temperature  ratio 

(12) 

0  =  T  /T^,  the  best  correlation  is  obtained  for  both  holes  and  slots, 
with  the  group  m”1'35  .  9  "°*68,  which  reduces  to  K  “1*35. 

The  best  overall  correlation  of  slot  and  holes  data  obtained  for 

a.  =  20°  (Figure  10)  is,  therefore,  obtained  against  the  group 
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Previous  film  cooling  data  for  values  of  a  =  10,  20,  30°  have  indicated  a 

2  _  i 

dependence  of  the  effectiveness  on  the  parameter  (Kq  cos  *  .  How¬ 

ever,  more  data  (over  a  wider  angular  range)  are  needed  to  validate  this 
correlation  fully. 

3.3.5.  Tests  and  Correlations  for  Slot-Injection  over  Curved  Surfaces 

The  test  channels  utilized  in  this  investigation  were  made  of 
plexiglass,  and  have  curbed  test  surfaces  scaled  to  the  following 
dimensions:  Test  Channel:  A  3”  convex  radius,  3  7/8"  concave  radius.  Test 
Channel  B:  7  " convex  radius,  7  7/8"  concave  radius.  Both  channels  have  a 
uniform  2"  x  7/8"  rectangular  cross-section  and  are  equipped  at  the 
downstream  end  with  a  choking  nozzle.  The  nozzle  has  a  1  1/2"  x  7/8" 
rectangular  throat,  which  produces  a  Mach  number  of  0.5  under  choked  flow 
conditions.  A  schematic  of  the  test  channel  is  shown  in  Figure  11.  The 
injection  slots  which  pass  the  coolant  from  the  plenum  chamber  to  the 
channel  are  2"  long  and  0.03"  wide.  The  angle  of  injection  is  20°  relative 
to  the  surface.  The  injection  slots  are  located  at  bend  angles  of  30°  in 
Channel  A,  and  17  in  Channel  B  from  the  channel  entrance.  These  locations 
assure  that  the  injection  is  made  in  region  of  nea  constant  static  pressure 
when  referred  to  Figure  12. 

Downstream  of  the  injection  slot,  each  curved  plate  is  instrumented 
with  14  thin  film  gauges:  8  along  the  centerline  to  measure  the  isothermal 
effectiveness,  and  6  across  the  plate  at  a  downstream  locations  x/s  =34. 

To  verify  the  flow  uniformity  a  map  of  the  distribution  of  the  gauges  is 
shown  in  Figure  13.  Each  channel  has  in  addition  a  Kistler  piezoelectric 
transducer  located  in  the  middle  of  the  side  wall.  Figure  14  gives  the 
isothermal  effectiveness  profile  on  the  3"  suction  surface  for  different 
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blowing  factors  which  shows  the  actual  increase  of  effectiveness  with 
increasing  m  and  the  usual  decrease  with  x.  Figure  15  summarizes  all  the 

=  x 

resu. ts  as  a  function  of  the  parameter  x  =  - ■=. 

sv 

These  results  show  clearly  that  a  smaller  radius  results  in  a  decrease 
of  effectiveness  on  the  concave  side  and  an  increase  on  the  convex  surface. 
The  values  of  the  effectiveness  on  the  convex  surface  exceed  that  on  the 
concave  surface. 

These  results  are  quite  expected  since  turbulent  mixing  spreads  across 

the  entire  boundary  layer  on  the  concave  side/1^  Further,  the  concave 

(13  14) 

side  promotes  the  formation  of  the  longitudinal  vortices  *  which 
increase  the  mixing  process. 

On  the  convex  surface,  the  thinner  boundary  layer  works  to  suppress 
radial  mixing.  The  shear  stress  appears  to  decrease  rapidly  from  the  near 
wall  region  to  the  edge  of  the  boundary  layer.  As  a  result,  the  turbulent 
heat  flux  tends  to  be  low  and  cooling  effectiveness  high.  When  the  effec¬ 
tiveness  data  are  plotted  against  the  correlating  parameter  A2  used  for 
flat  plate  data,  the  data  cluster  well  but  they  sho.»  scatter  particularly 
for  the  data  obtained  in  the  channel  with  a  short  radius  of  curvature. 

A  new  correlation  is  proposed  for  curved  walls.  It  includes  the 
effect  of  curvature  and  it  has  the  general  form 
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Y  =  K~2  [cos2a  ±  X  ~  sin2a]  | 

A  0  H  o 

Where  R  is  the  radius  of  curvature  of  the  wall,  X  is  an  empirical  non- 
dimensional  constant,  Dg  is  the  equivalent  slot  diameter  4A/P;  the  positive 
sign  applies  to  suction  (convex)  surface  data  and  the  negative  sign  to 
pressure  (concave)  surface  data.  Figure  16  shows  the  data  plotted  against 
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Y^  for  X  =  92.  The  data  cluster  well  but  a  small  divergence  exists  for  the 
data  obtained  on  the  concave  wall  of  the  channel  with  a  small  radius  of 
curvature . 


j.3.6  Steady  State  Experiments  on  a  Flat  Plate.  A  Validation  of  the 


Short  Duration  Technique 

The  correlations  developed  above  have  bean  obtained  in  short  duration 
tests  taking  great  care  to  obtain  steady  state  conditions  and  to  use  proper 
scaling.  In  the  scaling,  the  data  are  obtained  on  a  wall  operating  at  near 
ambient  temperatures  whereas  the  blade  surface  will  be  around  2000°R  in  an 
actual  turbine.  The  correlations  are  quite  reasonable  but  their  validity 
had  not  been  tested  at  the  conditions  of  an  actual  gas  turbine. 

The  validations  tests  were  conducted  in  a  steady  state  experiment 
using  a  flat  plate  cooled  by  a  double  row  of  holes  injecting  air  at  an 
angle  of  20°  to  the  combustion  gas  stream.  Test  conditions  were  varied  to 
provide  gas  temperatures  from  about  220O°R  to  3150°R.  This  corresponds  to 
gas-to-wall  temperature  ratios  of  3.6  to  5.2.  An  air  preheater  was  used  to 
heat  the  coolant  air  up  to  800°r.  This  corresponds  to  coolant  to  gas 
ratios  of  0.28  to  0.36  with  blowing  ratios  of  0.1  to  2.1.  Mach  numbers 
over  the  plate  were  kept  between  0.32  to  0.44.  The  heat  source  was  a 
combustor  burning  JP4  with  oxygen  diluted  witn  nitrogen  to  provide  desired 
gas  temperature.  The  combustor  and  the  facility  are  normally  used  for  MHD 
tests.  The  test  flat  plate  formed  one  of  the  four  walls  of  a  constant  area 
channel.  The  flat  plate  was  5"  wide  and  18"  long.  Upstream  of  the  flat 
plate  is  an  injection  block  with  the  two  rows  of  holes  having  a  diameter  of 
0.191".  The  pitch  between  holes  of  the  same  row  is  equal  to  two  diameters. 
Four  heat  transfer  gauges  made  of  three  thermocouples  located  at  different 


depths  are  located  along  the  longitudinal  axis  of  symmetry  of  the  flat 
plate.  The  gauges  are  located  at  distances  of  3*28",  5.53".  8.53"  and 
11.53"  downstream  of  the  second  line  of  holes,  corresponding  to  y/p  =  17.2, 
28.9,  44.7  and  60.4.  Figure  17  gives  a  schematic  of  the  film  cooling  test. 
It  indicates  the  location  of  pressure  sensors  and  thermocouples.  Figure  18 
gives  the  variation  of  isothermal  effectiveness  as  a  function  of  the 
blowing  ratio  at  the  different  axial  locations.  Figure  19  gives  the 
variation  of  measured  isothermal  effectiveness  with  the  correlation 
parameter  The  data  appear  to  be  as  well  correlated  as  they  were  for 

shock  tunnel  data.  Further,  the  comparison  (Figure  20)  of  all  shock  tunnel 
data  with  the  hot  long  duration  data  indicates  that  the  hot  data  fall 
within  the  correlation  hand  of  all  shock  tunnel  data.  This  good  agreement 
validates  the  simpler,  lower  cost  shock  tunnel  and  other  short  duration  low 
temperature  techniques. 
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SECTION  4 


4.  Analytical  Study  of  the  Thermo!  and  Fluid  Mechanical  Evolution  of  a 


Cooling  Film  Injected  from  a  Single  Line  of  Inclined  Round  Holes 


4 . 1  Introduction 

The  use  of  a  film  of  cooling  air  is  a  promising  method  of  reducing  the 
heat  load  and  metal  temperature  of  hot  gas  turbine  components  exposed  to 
very  high  temperature  combustion  gases. 

The  rate  of  mixing  of  coolant  with  the  hot  main  stream  determines  the 
effectiveness  of  the  film  as  a  thermal  barrier.  The  film  of  coolant  air 
can  be  obtained  by  blowing  air  through  a  porous  surface,  a  slot,  or  an 
array  of  holes. 

Transpiration  cooling  uses  relatively  small  cooling  flows.  The  momen¬ 
tum  of  blown  air  is  negligible  compared  to*  the  momentum  of  the  main  flow 
and  boundary  layer  theory  can  be  used  to  determine  the  thermal  effective¬ 
ness  of  the  film.  However,  the  technological  implementation  of  trans¬ 
piration  cooling  has  met  severe  difficulties  and  in  practice  the  film  is 
obtained  by  injection  slots  and  more  often  through  arrays  of  holes. 

In  a  few  instances,  the  film  can  be  injected  through  slots  parallel  to 

the  surface.  In  these  cases,  the  mixed  layer  is  dominated  by  wall  shear 

(15) 

stress  and  the  evolution  of  the  film  can  be  described  analytically 
(16) 

This  analysis  considers  a  more  general  type  of  film  injection  using  a 
line  of  round  holes  inclined  at  a  certain  angle  to  the  surface.  (The 
projection  of  the  hole  axis  on  the  surface  is  oriented  in  the  streamwise 
direction).  In  this  case,  boundary  layer  theories  cannot  be  applied,  since 
the  momentum  normal  to  the  wall  is  not  negligible  compared  to  tne  main 
stream  momentum.  Further,  the  jets  separate  from  the  wall  at  the  point  of 
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injection  and  a  pair  of  contra-rotating  vortices  is  created  within  each  jet 

by  the  velocity  field  interaction.  The  strength  of  the  vortices  can  be 

determined  from  this  interaction.  The  achievement  of  good  film  cooling 

effectiveness  requires  the  reattachment  of  the  individual  jets  to  the  wall 

at  a  distance  which  is  short  compared  to  the  length  of  surface  which  has  to 

be  effectively  cooled.  The  latter  consideration  puts  a  limit  on  the  normal 

momentum  of  the  jets.  The  decay  of  these  contrarotating  vortices  is  found 

to  be  negligible  over  the  length  of  wall  to  be  effectively  cooled.  Since 

the  holes  are  fairly  closely  spaced,  there  is  a  mutual  interaction  between 

vortices  issued  from  neighboring  holes.  The  mixing  of  the  two  streams  is 

due  to  both  turbulent  diffusion  and  vortical  entrainment. 

Experimental  work  has  shown  the  importance  of  the  separated  region 

downstream  of  the  separation  point  and  of  the  contrarotating  vortices  in 
( 17) 

the  mixing  process  .  A  serai-empirical  analysis  of  mixing  by  turbulent 

diffusion  and  vortical  entrainment  was  developed  for  the  case  of  inclined 
( 18 ) 

slot  injection  .  The  concept  of  turbulent  and  vortical  mixing  for  films 

injected  from  a  line  of  inclined  round  holes  was  first  used  in  an 

analysis^)  which  determined  the  trajectory  of  a  single  jet  and  limited 

the  effect  of  vortical  entrainment  to  a  short  distance  (15  hole  diameters) 

compared  to  length  of  effectively  cooled  surface. 

(20) 

The  present  analysis  applies  this  mixing  process  and  vortex  inter¬ 
action  to  the  initial  trajectory  of  the  jets  up  to  the  formation  of  a 
continuous  layer  and  to  the  growth  of  the  mixed  layer. 

The  analysis  is  divided  into  two  parts.  The  first  describes  the 
trajectory  and  mixing  of  the  jets  and  their  reattachment,  whereas  the 
second  determines  the  growth  of  the  mixed  layer  and  the  change  of  crosswise 
average  film  cooling  effectiveness  with  distance  from  the  injection  point. 
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A  comparison  of  the  present  theory  is  made  with  available  data  of  film 
effectiveness  measured  for  films  injected  from  a  single  line  of  inclined 
holes. 


4.2  Analysis  of  the  Jet  Trajectory  and  Mixing 


Jets  issue  from  a  row  of  round  holes  (spacing  s  and  diameter  D^)  with 


an  average  velocity  u^q  and  a  density  p ^  at  an  angle  Bq  (Fig.  21)  and 


penetrate  into  a  main-flow  of  velocity  u^  and  density  p  outside  the 
boundary  layer.  The  flow  field  is  steady. 

The  boundary  layer  of  the  main  flow  is  assumed  to  be  turbulent  with  a 
thickness  6: 


1/7 


,Z. 

u  =  (-)  u 

if 

z  < 

6 

0  CO 

• 

U  s  U 

if 

z  > 

<5 

(1) 


Contrarotative  vortices  appear  inside  the  jets  and  are  the  characteristic 
feature  of  the  flow  field.  A  qualitative  description  of  these  vortices  was 
given  in  Reference  21  and  the  first  quantitative  description  of  them  is  in 
Reference  22. 


4.2.1.  Equations  of  Motion 

The  equations  are  written  in  curvilinear  coordinates  (Fig.  24).  Each 


cross  section  of  the  jet  has  an  area  A.,  an  average  velocity  u.,  an  average 

J  J 


density  p.  and  an  average  temperature  T^;  8  is  the  angle  of  the  jet  with 

J  J 


the  horizontal  (Fig.  22);  E  is  the  entrainment  of  the  main  flow  per  unit 
length  of  the  jet.  u,p  and  u^  are  the  tangential  and  normal  component  of 
the  velocity  of  the  main-flow  relative  to  the  jet  in  the  vicinity  of  the 
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jet  (main-flow  velocity  induced  by  the  vortices),  R  is  the  local  radius  of 
curvature  of  the  jet,  cD  is  the  drag  coefficient  of  the  jet,  and  D  is  the 


large  diameter  of  the  jet. 


The  equations  of  motion  are: 


.(17) 


Equation  of  Continuity 


dA  (pjUjAj}  £  E 


Momentum  Equation 


h  4  *  eut 


Force  Equation  (normal  to  the  jet) 


V  1  2 

3  Arr  =  EUK  *  2  D 


(The  centrifugal  force  balances  the  normal  components  of  the  drag  acting  on 
the  jet  and  of  the  variation  of  momentum  due  to  the  entrainment.) 

The  force  equation  can  also  be  written: 

-»jVj  *  Eu»  *  1  »cduk  D 

A  dilution  parameter  in  the  jet,  Uj,  is  defined  as  the  ratio  of  main  flow 
which  has  been  entrained  inside  the  jet  to  the  initial  mass  flow  of  the 
jet.  Therefore: 


dA  "  p-  U.  A. 

PJO  jo  jo 


The  average  temperature  in  the  jet  is  given  by  the  equation  of  energy: 


T.  +  ucT 

T  J°  S 

j"  1  +  M. 


If  the  flow  is  assumed  incompressible, 
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Equations  (2),  (3),  (5),  (6)  and  (8)  are  solved  in  a  step-by-step  calcu¬ 
lation. 


4.2.2  Entrainment 

The  entrainment  of  the  main-flow  is  due  partly  to  the  contrarotative 
vortices  generated  inside  the  jets  and  partly  to  turbulent  diffusion. 

According  to  Ref.  19  the  entrainment  due  to  the  vortices  is  equivalent 
to  a  sink  singularity  concentrated  just  behind  the  jet. 

The  complex  potential  of  a  flow  due  to  a  sink  of  strength  E  located  at 
z  =  a  and  with  a  velocity  u^  at  infinity  is: 


f .  (z)  =  UMz  -  ~ —  L 
l  N  2irp 


z  =  x  +  iy 

If  the  cross  section  of  the  jet  remains  almost  circular,  the  flow  has  to  be 
such  that  the  circle  |z|  =  a  is  a  streamline  (Fig.  23).  Hence  the  complex 
potential  of  the  flow  in  presence  of  a  circular  cylinder  of  radius  a  is 
(Ref.  24,  p.  422:  Image  effect  of  the  sink  relative  to  the  circle): 

2  2 

f(z)  =  U  (z  ♦  — )  -  Log  [~Z~~]  +  const.  (10) 

z  2irp  Z 

The  force  acting  on  a  unit  length  of  the  jet  is  given  by  Blasius 
(5) 

formula  : 
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(11) 


dF 

ds 


<||)  dz  =  EUn 


This  force  can  also  be  expressed  in  terms  of  a  drag  coefficient  c^: 


ds  ‘  2  CDD 


From  Eqs.  (11)  and  (12) 


E  =  r  V 


E  =  pE^D 


with  E^  = 


CD 

2 


(12) 


(13) 


This  relation  between  the  coefficient  of  entrainment  E1  and  the  drag 

coefficient  cD  is  assumed  to  be  valid  even  when  the  cross  section  of  the 

jet  does  not  remain  circular.  In  fact  the  cross  section  of  the  jet  is  . 

changed  from  the  circle  to  a  characteristic  kidney  shape  :  the  sides  of 

the  jet  are  shifted  downwards  faster  than  the  core  of  the  jet  because  of  a 

larger  mixing  with  the  main  flow  (Fig.  24).  In  the  calculation  it  is 

assumed  that  the  cross  section  has  an  elliptic  shape.  The  drag  coefficient 

(O'3,) 

of  an  elliptic  cylinder  is  about  1.8  •  Hence  E1  =  c^/2  =  0.9. 

The  entrainment  due  to  turbulent  diffusion  is  proportional  to  the 

difference  of  tangential  velocities  and  to  the  circumference  of  the  jet  C: 

E  =  PE2|uT-uj|C  (14) 

( 26 ) 

According  to  Ref.  25,  E2  =  0.08.  Similar  measurements  lead  to 
E2  =  0.116.  It  is  assumed  that  E^  =  0.08. 

Therefore  the  entrainment  of  the  main-flow  is 

E  =  pE1  u^D  +  pE2|u^,  -  u j  ( C  (15) 

with  E1  =  0.9  and  E2  =  0.08 
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4.2.3  First  Section  of  the  Jets 


When  the  jets  exhaust  from  the  holes,  the  lower  surface  of  the  jets  is 
touching  the  wall  (Fig.  25).  In  this  section  the  main  velocity  is  assumed 
to  be  unchanged,  therefore: 
u-  =  u  cos  8 

u^  =  u  sin  8  (16) 

In  fact,  there  is  a  contraction  of  the  main  flow  due  to  the  jets,  and  the 
main  velocity  is  increased  slightly.  This  contraction  was  neglected  in  the 
calculation.  In  section  1  the  jet  is  turned  from  bq  to  8^  and  is  squeezed 
between  the  main  flow  and  the  wall  so  that  the  shape  of  the  cross  section 
becomes  elliptical.  The  calculation  of  the  shape  of  the  cross  section  is 
carried  out^20\ 

4.2.4  Creation  and  Strength  of  the  Contrarotative  V  ^rtices 
In  Section  4.2.2  of  the  jet  (Fig.  25),  the  jet  is  unattached.  The 
component  of  the  main-flow  velocity  normal  to  the  jet  in  cross  section  AB 
of  Fig.  25  is  Uj|  :  u  sin  B 1 .  Two  contrarotative  vortices  appear  inside  the 
jet  to  balance  this  normal  velocity  so  that  the  boundary  of  the  jet  is  a 
streamline  (Fig.  26).  B  and  D  arv,  the  small  and  large  diameter  of  the 
elliptical  cross  section  of  the  jet.  The  stream  function  in  the  first  free 
jet  cross  section  is  such  that 


V 


3f_ 
3  z 


and  U 

o 


dj/_ 

3  r 


Therefore 


r  r  r  r 

»  ■  »y  ■  k  to. rW  -vn-' 

2*  r2rnrg.  r2n 


(17) 
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where  it  is  assumed  that  the  velocity  induced  by  a  curved  vortex  line  is 
the  same  as  the  one  induced  by  an  infinite  rectilinear  vortex  line  tangent 
to  the  curved  vortex  line  (i.e.  the  flow  starts  rotating  inside  the  holes). 
By  symmetry,  V ^  =  0 

The  boundary  of  the  jet  is  a  streamline,  thus  =  0: 


From  Eq.  (19): 


s  -  Y„ 


S  +  Y- 


r  =  tfU^/[  2  2 

N  yt  +  r 

1  4 


p  2  2  2 

(s-Y.)  +  B  (s+YT)  +  B 
1  4  1  4 


ns  -  Y„ 


ns  +  Y_ 


(ns-YT)2+  B2  ( ns+YT)2+  B2 
1  4  1  4 


The  values -of  r  obtained  from  this  calculation  are  close  to  the  ones 
obtained  from  the  empirical  relation  given  in  Reference  22: 

r  =  1.4  u.0  sin  bq  Dq  (2 

where  it  is  assumed  that  only  the  normal  component  u^q  sin  Bq  of  the  jet 
velocity  induces  some  vorticity. 

In  the  configuration  studied  in  Ref.  27,  the  calculation  leads  to  an  a 


dimensional  strength  of  the  vortices: 


=  .490 


while  the  application  of  Ref.  22  empirical  relation  gives: 


=  .474 


The  difference  is  only  3%. 


4.2.5  Induced  Normal  Velocity 

Reference  20  shows  that  the  damping  of  the  vortices  due  to  viscosity 
is  felt  from  a  distance  of  about  400  hole  diameters.  Therefore,  in  the 
range  usually  considered  in  film  cooling  (up  to  100  hole  diameters)  the 
damping  can  be  neglected;  the  fluid  is  assumed  inviscid  except  for  the 
entrainment  of  main-flow. 


According  to  Kelvin’s  circulation  theorem,  the  circulation  around  any 
closed  material  curve  is  invariant  in  an  inviscid  fluid.  This  can  be 
applied  to  curve  ;  in  Fig.  27  which  encircles  half  of  the  jet  cross 
section. 

By  Stoke' s  theorem,  the  vorticity  included  inside  the  curve  c  is 
invariant;  the  strength  of  the  contrarotative  vortices  is  invariant. 

Point  A  at  the  top  of  the  jet  cross  section  is  a  stagnation  point  in 
the  cross  section  (Fig.  28).  With  the  notations  of  Fig.  28; 


UA  = 


„  r  r _ _T 

UN+  IT  v  2 


s  -  Y„ 


s  +  Y_ 


Y"  +  (Bc-  2„) 
1  2  1 


(s-Yt)+(B2-  zJ  ( s+Y~ ) B  -  ZT)2 

2  1  2 


ns  -  Y„ 


ns  +  Y„ 


(ns-YT)2+  (B  -  ZT)2  (ns+Y2)  +  (B  -  ZT)2 


-]  =  0 


Therefore,  the  normal  velocity  of  the  main  flow  induced  by  the  vortices  in 


the  vicinity  of  the  jets  is: 


ns  +  Y_ 


ns  -  Y„ 


U  -  [ - 1 - 

N  *  Yt2+(B-zt) 


+  z  (- 


2  m  (s-Yt)+(B  -  Z— ) 2 
1  2  1 


(s+Yt)2(B  -  ZT)2 
2 


-)]  (23) 


From  the  conservation  of  mass  flow,  the  horizontal  velocity  of  the  main 
flow  is  unchanged,  neglecting  the  increase  of  the  jet  cross  section  area. 
Consequently,  the  velocity  of  the  main  flow  is  deflected  downwards  from  u 


to  u'  such  that  the  normal  corapont  is  u^  (Fig.  29).  The  tangential 


velocity  is 


= 


T  "  cos  6  N 


-  U„tan8 


(24) 
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4.2.6  Trajectory  and  Growth  of  the  Jets 

The  jet  centerline  is  obtained  from  the  equations  of  motion  (§  1)  with 
the  entrainment  given  in  (§  2)  and  the  normal  and  tangential  components  of 
the  main  velocity  given  in  (§  5).  Due  to  the  normal  velocity  induced  by 
the  vortices,  the  jets  are  deflected  downwards  below  the  horizontal.  The 
jets  reattach  to  the  wall  along  line  R  (Fig.  30).  The  shape  of  the  jets 

remains  the  same  as  long  as  they  are  free  jets.  If  the  jets  reattach  to 

the  wall  before  touching  one  another,  it  is  assumed  that  they  continue  to 
grow  only  in  the  lateral  direction  (Fig.  31a).  If  the  jets  are  in  contact 

with  one  another  before  reattaching,  it  is  assumed  that  they  continue  to 

grow  only  in  the  vertical  direction  (Fig.  31b). 


4.2.7  Trajectory  of  the  Vortices 

In  a  cross  section  of  the  jet  (Fig.  32)  the  vertical  velocity  of  a 
vortex  located  at  (Y^,  ZT)  relative  to  the  jet  centerline  in  curvilinear 
coordinates  (Fig.  24)  is: 


V  ,,  X  1  1  1 

vz  "  "  n  +  2*  2Y  "  S-2Yt  +  S+2YT 


nS-2YT  +  nS+2YT 


For  a  step  As  along  the  jet  centerline,  the  step  in  time  is  At  = 


As/Uj  and 


AZ  =  Y  AT  =  -—  AS 

j 


From  the  position  of  the  vortices  at  a  given  step,  Eq.  (23)  gives  the  value 
of  and  Eqs .  (25)  and  (26)  lead  to  the  position  of  the  vortices  at  the 
next  step.  Since  the  location  of  the  vortices  in  the  first  free  jet  cross 
section  is  known  [Z^  =  0  and  computed  from  Eq.  (20)],  the  path  of  the 
vortices  is  deduced  from  the  trajectory  of  the  jet  centerline. 
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4.2.8  Cooling  Effectiveness  for  Reattached  Jets 

When  the  jets  are  touching  the  wall,  it  is  assumed  that  the  laterally 
averaged  wall  temperature  is  the  same  as  the  averago  temperature  of  a  layer 
with  a  thickness  equal  to  the  small  diameter  of  the  jets.  When  the  jets 
are  in  contact  with  one  another  and  are  touching  the  wall,  they  form  a 
layer  over  the  plate  (Fig.  33). 

The  dilution  parameter  at  the  beginning  of  the  layer  is 

A-A 

"jF  *  “j  *  —  *“j>  W 


The  adiabatic  effectiveness  of  the  cooling  is  defined  as 


n  =  t-T 

jo 


with 

T  s  x  s  TJ°  * 

w  Lo  1  +  Ujp 


therefore 

1  _^j  1 

n  ’  ^JF‘  A  1+»j 
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4 . 3  Analysis  of  the  Mixed  Layer 

When  the  jets  are  in  contact  with  the  wall  and  with  one  another,  the 
cool  flow  constitutes  a  layer  over  the  plate.  For  the  calculation  it  is 
assumed  that  the  initial  thickness  of  the  layer  is  equal  to  the  small 
diameter  of  the  last  cross  section  of  the  jets. 


4.3.1  Equations  of  the  Layer 


An  additional  dilution  of  the  main  flow  in  the  cool  flow  occurs  in  the 
layer.  A  dilution  parameter  in  the  layer  y^  is  defined  as  the  ratio  of 
main  flow  which  has  been  entrained  inside  the  layer  to  the  initial  mass 
flow  of  the  layer. 

Equation  of  continuity: 


PLULh  =  (1  °I.ouLoho 


(30) 


Conservation  of  momentum: 


0  *  »L>  UL  =  UL0  *  »L  U-  °r  = 


u_  +  u,  U 
Lo  « 


L  L  1  +  u 


(3D 


Equation  of  energy: 


<>V  V  Lo*  »lt 


or 


T  = 


T,  +  y .  T 
Lo  L 

1  +  y. 


(32) 


If  the  flow  is  assumed  incompressible,  the  equation  of  energy  can  be 
written : 


,  1  1  L 

( 1  +y .  )  —  =  -  +  — 

L  p  L  pLc  p 


(33) 


or 


pLop 

p.  =  ( 1  +y .  )  - - — 

L  L  p  +  y .  p 


LK  Lo 


Substituting  Eqs.  (31)  and  (33)  into  the  equation  of  continuity  (30) 
leads  to: 

p. 


1  +  m 


Lo 


h  =  ( 1  +y  L ) 


l  P 


1  +  u 


U 

a 

L  TT 


(34) 
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4.3.2  Dilution  Due  to  the  Vortices 


The  flow  is  considered  inviscid.  The  velocity  at  the  wall  has  to  be 


tangent  to  the  wall.  Therefore,  the  velocity  field  is  the  same  as  the  one 


induced  by  the  vortices  and  their  images  symmetrical  to  the  wall  (Fig.  31*)  • 


The  mass  flow  of  main  fluid  entering  the  layer  per  unit  area  due  to 


the  vortices  only  is  the  same  as  the  mass  flow  which  would  leave  the  layer 


in  the  velocity  field  induced  by  the  vortices  and  their  wall  images: 


♦  =  2  -  ;  V  (M)  dy 
o 


where  a  is  the  point  of  horizontal  velocity  at  the  boundary  of  the  layer, 

.  .  ...  _ , _ „  (20) 


With  the  notations  of  Fig.  14,  and  assuming  a  *  Y^ 


of  r,..  r  .  „  r((n+1)s)  +  (  r) 

—  Llog  r—  +  l  log  [ - = - r—  X 

irS  n  2r  n=0  ((n+1)s)  +(h-z) 

P 


[(ns+2YT)2+  (h-zT)]2[((n+1)s-2YT)2  +  (h-zr)2] 


[(ns+2YT)2  +  (h+zT)2]  C((n+1)s-2YT)2+  (h+zT>2] 


4.3.3  Dilution  Due  to  Turbulence 


Acording  to  Fick's  law,  the  mass  flux  entrained  into  the  layer  by 


turbulence  is: 


,t  =  -  P°t  8rad  c  dx 


where  c  is  the  concentration  of  main  fluid,  grad  c  =  -ac/az 


and  Dfc  is  the  turbulent  diffusion  coefficient, 


h  „ 
p  T—  U 


K(n)  dx 
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where  h^  is  the  distance  between  the  wall  and  the  inflection  point  of  the 
concentration  profile,  4  is  the  overall  dilution  parameter  defined  as  the 
ratio  of  fluid  originating  from  the  main  flow  to  the  fluid  originating  from 
the  holes: 

"  =  "jF  *  UL  <’•  *  ujF>  *  “jF  +  "L  *■  Vl  (39) 

and 

K(n)  =  0.052353  -  0.032051n  +  0.006582n2  -  0.000433n3  (40) 

Overall  Dilution 

In  a  step  ax  the  mass  flow  of  main  fluid  entering  the  layer  by 
vortical  and  turbulent  mixing  is 

Aid 

Aij»e=  4  +  ~ ax  (41) 

where  $  is  given  by  Eq.  (36)  and 
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dm 
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Ax 


dx 


by  Eqs .  (38)  and  (40) . 


The  variation  of  the  layer  dilution  parameter  is  such  that 
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(42) 


4.3.5  Velocity  Profile 


The  average  velocity  in  the  layer  u^  is  given  by  the  conservation  of 
momentum  (31)  neglecting  the  wall  friction.  The  effect  of  the  wall 
friction  is  to  induce  a  velocity  profile  which  is  assumed  to  be  of  the  form 


U  =  U  (£) 
00  n 


1/n 


(43) 
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When  the  mixed  layer  is  a  boundary  layer  of  thickness  h,  the  velocity 
profile  is  assumed  to  be  a  turbulent  one: 


U  =  U  (-£)  (44) 

oo  n 

It  is  assumed  that  the  ratio  of  average  velocities  deduced  from  Eqs.  (43) 


and  (44)  is  equal  to  u./u  : 

L*  ® 


'L  l  <l>  dz  l  r  ”dt  8 

h  1  s  V  1  =  1. 

/  <f>  fn  J  r  7dt  7(1 


Therefore 


UL/U- 

8/7-00“ 

L  « 


Substituting  the  expression  (31)  of  u^  into  (45): 


7^%t 

V  8~7  ULC 


If  u^q  <  u^,  n  is  increasing 


from  r.0  = 


8-7 


when  u  is  ve^y  large  (Fig.  36),  and  if  u  >  u  n  is  decreasing  from  n  to 

L  LU  »  U 


4.3.6  Trajectory  of  the  Vortices 

The  fluid  inside  the  layer  cannot  leave  the  layer.  Hence,  the 
vertical  velocity  of  a  fluid  particle  inside  the  layer  near  the  upper 


surface  of  the  layer  cannot  be  larger  than  the  growth  rate  of  the  layer 
which  is  very  small.  Therefore  the  velocity  of  this  fluid  particle  in  the 
cross  section  is  almost  horizontal. 

In  order  to  have  an  horizontal  velocity  at  the  surface,  the  velocity 
field  inside  the  layer  is  the  same  as  the  one  induced  by  the  vortices  and 
their  symmetric  images  relative  to  the  wall  and  relative  to  the  upper  line 
of  the  layer,  and  also  by  the  symmetric  images  of  those.  The  velocity 
field  is  induced  by  an  infinity  of  double  rows  of  vortices  at  a  distance  2h 
from  one  another  (Fig.  37).  The  velocity  of  a  vortex  located  at  (YT>  Z T> 
due  to  the  other  vortices  is  derived. 

Knowing  the  induced  velocity,  the  vortex  trajectory  can  be  computed  as 
shown  in  Fig.  38.  Since  the  layer  is  growing,  the  vortex  is  on  an 
ascending  spiral  inside  the  layer  (Fig.  39). 

14.3.7  Film  Cooling  Effectiveness 

The  adiabatic  film  cooling  effectiveness  of  the  layer  is  defined  as 


where  Hq  is  the  effectiveness  at  the  beginning  of  the  layer.  Because  of 
the  presence  of  the  wall,  the  diffusion  of  fluid  particles  is  not 
symmetrical  and  there  is  a  reflection  .  the  wall  (Fig.  *40) .  Therefore, 
the  concentration  and  temperature  profiles  have  an  inflection  point  at  a 
distance  h-j.  from  the  wall  which  is  moving  upwards. 
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When  h  is  very  large  hj  *  h/2.  At  the  beginning  of  the  layer  it  is 
assumed  that  hj  *  0.7  hQ.  Therefore, 

hj  *  0.2  hQ  +  0.5  h  (48) 

The  temperature  profiles  are  shown  in  Fig.  41.  The  average 
temperature  in  the  layer  is  approximated  by 


T  + 


w 


(1 


+ 


(49) 


and  the  adiabatic  effectiveness  is: 
u  T  -  T, 


n  = 


hT  T  -  T,  o 
I  Lo 


replacing  T^  by  its  expression  (32) 


h  1 

"  S  hi  1  *uL  n° 


(50) 


This  analysis  can  be  extended  to  multiple  lines  of  equally  spaced 
holes  using  the  additive  superposition  proposed  in  Ref.  29.  The  validity 
and  limitations  of  this  treatment  are  determined  in  Refs.  30  and  31. 


4.4  Comparison  with  Experimental  Data 

4.4.1  Jet  Trajectory,  Growth  of  the  Layer  and  Path  of  the  Vortices 
The  calculation  has  been  made  for  the  cases  presented  recently  in 
Reference  27.  It  includes  the  results  presented  in  Ref.  32.  The  spacing 
between  the  holes  is  3  hole  diameters,  the  angle  of  injection  is  35  degrees 
and  the  density  velocity  ratio  defined  as 


m  = 


p  .  U. 
jo  jo 


is  0,5. 


(51) 
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Since  the  density  ratio  p  jQ/p  is  about  0.85,  the  velocity  ratio 

U.-/U  is  about  0.59. 
ju  ® 

When  the  ratio  of  the  boundary  layer  thickness  to  the  hole  diameter  is 
0.245,  the  jet  trajectory  and  the  path  of  the  vortices  inside  the  jets  are 
presented  in  Fig.  42.  Th  reattachment  to  the  wall  is  found  at  a  distance 
of  4.35  hole  diameters  from  the  holes.  The  growth  of  the  layer  and  the 
path  of  the  vortices  inside  the  layer  are  presented  in  Fig.  43  for  the  same 
case.  At  the  bottom  of  the  figure  the  projection  of  the  path  of  the 
vortices  in  the  cross  section  of  the  layer  shows  that  the  vortices  have  an 
helicoidal  motion  inside  the  layer. 

4.4.2  Film  Cooling  Effectiveness 

The  calculated  film  cooling  effectiveness  compared  to  the  experimental 
results  for  the  case  described  above  is  presented  in  Fig.  44.  The 
calculated  curve  is  above  the  experimental  results  of  average  film  cooling 
effectiveness  near  the  reattachment,  but  the  agreement  is  good  at  a  certain 
distance  from  the  reattachment. 

The  calculation  has  also  been  done  for  some  cases  presented  in  Hof. 

33.  The  geometric  configuration  is  the  same  as  for  the  case  of 

Ref.  27  (s  r  3Dq  and  S  =  35°).  Two  curves  are  presented  in  Fig.  45,  where 

the  ratio  of  the  boundary  layer  thickness  to  the  hole  diameter  is  0.174: 

one  for  the  density  velocity  ratio  m  =  0.31  (p  .A/p  =  0.742  and  u.n/u  =  0.418) 

ju  ju 

and  the  other  for  m  =  0.70  (p  .n/p  =  0.769  and  u.A/u  =  0.910).  In  these  cases 

ju  ju  « 

the  agreement  between  the  calculation  and  the  experimental  results  is  good 

except  very  near  the  injection  (up  to  5  hole  diameters). 

( ■34 ) 

Two  cases  presented  recently  for  a  single  line  of  holes  were  also 
calculated.  The  spacing  between  the  holes  is  3  hole  diameters,  the 


ejection  angle  is  45°,  the  ratio  of  the  boundary  layer  thickness  to  the 
hole  diameter  is  0.104,  and  the  density  ratio  is  p jQ/p  =  0.94.  Two  curves 
are  shown  in  Fig.  46:  one  for  a  density  velocity  ratio,  m  =  0.2T 
(u^/u  =  0.29),  and  the  other  for  m  =  0.54  (u^/u^  =  0.57).  The  agreement 
between  the  calculation  and  the  experimental  results  is  good  except  very 
near  the  injection  point.  The  model  is  not  very  accurate  near  the  in¬ 
jection  because  the  velocity  profile  inside  the  tubes  has  been  assumed 
constant.  In  fact,  the  low  momentum  fluid  is  deflected  much  more 
rapidly  than  the  core  fluid  by  the  main  flow  and  increases  the  cooling 
effectiveness  near  the  injection  point. 

4.4.3  Distance  of  Reattachment 

The  calculated  distances  xR  between  the  holes  and  the  reattachment 
point  in  the  different  cases  presented  above  are  as  follow: 

Ref.  27 

§-  =  3.  Bo  =  35  m*  0-5 

o 


m  =  0.70 


R 


=  12.70 
o 


Ref.  34 


=  3. 


6o  =  45 


-g-  =  0.104 
o 


ra  =  0.27 


=  2.20 


m  =  0.54 


kR 


=  3.86 


These  results  show  that  the  distance  between  the  holes  and  the  reattachment 
point  increases  when  m  increases  (more  penetration  of  the  jets)  and  also 


when  6 /DQ  increases  (less  deflection  of  the  jets). 


This  study  provides  a  comprehensive  analysis  of  the  fluid  mechanical 
and  thermal  evolution  of  the  film  formed  by  the  coolant  injected  through  a 
line  of  closely  spaced  round  holes  inclined  to  the  surface. 

The  first  part  of  the  analysis  describes  the  growth  of  the  jets  from 
the  injection  to  the  formation  of  a  continuous  mixed  layer.  The  strength 
of  the  vortices,  the  mutual  interaction  of  vorcices  induced  in  neighboring 
jets,  and  the  mixing  due  to  Mth  turbulent  diffusion  and  vortex  entrainment 
are  used  to  determine  the  trajectory  and  growth  of  the  jets.  The  second 
part  of  the  analysis  describes  the  growth  of  the  mixed  layer  as  influenced 
by  both  turbulent  diffusion  and  vortex  entrainment  since  the  decay  of  the 
vortex  strength  is  found  to  be  negligible  over  the  length  of  wall  effec¬ 
tively  cooled. 

The  results  of  the  analysis  compare  very  satisfactorily  with  available 
experimental  data  of  average  film  cooling  effectiveness  obtained  for  a 
single  line  of  holes. 
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The  analysis  can  be  extended  to  the  film  produced  by  Injection  of 
coolant  through  multiple  lines  of  holes  using  the  principle  of  super- 

(pQ  ) 

position  within  the  limitations  determined  in  References  30  and  31. 
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4 . 5  Nomenclature 


w 


E1’E2 


f 

* 


jF 


U 

n 

r 


m 


e,t 


m 


m,c 


Q 


m 


Temperature  in  the  layer  Cayerage) 
Wall  temperature  (adiabatic) 
Entrainment  of  the  main  flow 
Entrainment  coefficients 


Drag  coefficient 

Complex  potential 

Stream  function 

Dilution  parameter  in  the  jet 

Dilution  parameter  at  the  beginning  of  the  layer 
Dilution  parameter  in  the  layer 
Overall  dilution  parameter 
Adiabatic  film  cooling  efficiency 
Strength  of  the  vortices 
Mass  flow  entrained  by  the  vortices 
Mass  flow  entrained  by  turbulent 
Turbulent  diffusion  coefficient 
Concentration 
Turbulent  viscosity 
Velocity  mixing  length 
’Concentration  mixing  length 
Turbulent  Schmidt  number 
Heat  transfer 
Heat  transfer  coefficient 


Density  velocity  ratio  m  ** 


Pi0U10 

PU 
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Nomenclature  (Cont'd) 


Subscripts: 


&  1 


Superscript: 


j  in  the  jets 

jO  at  the  beginning  of  the  jets 

L  in  the  layer 

LO  at  the  beginning  of  the  layer 

v  at  the  vail 

T  relative  to  the  vortices 
iso  isothermal 
a  adiabatic 

c  coolant 

R  rcattichment 

average 


SECTION  5 


5.  Aerodynamic  performance  and  Heat  Transfer  Characteristics  of  Film- 
Cooled  Airfoils  in  a  Subsonic  Cascade  Flow 

5 . 1  Introduction 

The  use  of  a  film  of  air  coolant  is  an  effective  technique  for  the 
reduction  of  heat  transfer  to  hot  components  of  gas  turbines.  Therefore, 
the  interaction  and  mixing  of  the  coolant  with  the  high  temperature  main¬ 
stream  over  the  blade  surface  influences  not  only  the  thermal  effectiveness 
of  the  film,  but  also  the  aerodynamic  profile  losses  of  film 
cooled  blades. 

The  film  of  coolant  can  be  obtained  by  blowing  air  through  a  porous 
surface,  arrays  of  holes,  and  slots.  (The  last  technique  is  seldom  used  in 
blade  cooling  because  of  mechanical  difficulties.)  The  first  technique, 
transpiration  cooling,  uses  cooling  flows  of  low  momentum  which  can  be 
described  by  boundary  layer  theory.  The  aerodynamic  and  heat  transfer 

Mr  £  ■37} 

performance  of  transpiration  cooled  blades  has  been  investigated v  0  ’ 

The  difficulties  of  technological  implementation  of  transpiration  cooling 
and  the  significant  aerodynamic  penalty  have  impeded  the  application  of 
this  mode  of  cooling. 

In  present  practice,  the  film  is  obtained  by  the  injection  of  a  film 
through  an  array  of  holes  inclined  to  the  wall.  The  cool  jets  issued  from 
these  holes  were  found  to  have  a  momentum  normal  to  the  wall  which  was  not 
negligible  compared  to  the  momentum  of  the  mainstream,  and  boundary  layer 
could  not  be  applied  to  the  development  of  the  film.  Further,  the  jets 
separated  at  the  point  of  injection  and  a  pair  of  contrarotating  vortices 
was  created  within  each  jet  by  the  velocity  field  interaction.  When  the 
jets  reattached  to  the  wall  and  formed  a  mixing  layer,  the  vortices 
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remained  imbedded  in  the  layer.  The  evolution  of  the  coolant  film  was, 
therefore,  determined  by  entrainment  of  the  main  flow  due  to  two 
mechanisms:  turbulent  diffusion  and  the  velocity  field  induced  by  the 
vortices.  The  concept  of  turbulent  and  vortical  mixing  for  cooling  films 
has  only 

( '18  8Q ) 

recently  been  formulated  for  the  evolution  of  films  over  flat  plates  ’ 
Although  considerable  information  exists  on  the  thermal  effectiveness 
of  films  issued  from  arrays  of  holes  over  flat  surfaces,  very  limited 
information  exists  on  the  effects  of  flow  curvature,  pressure  gradients  and 
turbulence  levels  on  the  behavior  of  a  coolant  film  along  a  turbine  blade. 
The  thermal  performance  of  a  blade  cooled  through  a  discrete  number  of 
lines  of  round  holes  has  been  investigated and,  more  recently,  the 
aerodynamic  and  heat  transfer  performance  of  a  full  coverage  film  cooled 

(41 ) 

blade  has  been  reported 

This  study  presents  a  comprehensive  set  of  aerodynamic  and  heat 

transfer  studies  obtained  on  film  cooled  blades  with  straight  suction  back 

(42) 

in  a  linear  cascade.  The  influence  of  location  of  the  coolant 
injection  on  performance  is  investigated. 

5.2  Experimental  Conditions  and  Apparatus 

The  aerodynamic  and  heat  transfer  measurements  reported  in  this  paper 
were  made  on  the  M.I.T.  hot  turbine  blowdown  facility  (Figure  47).  This 
facility  uses  a  pebble  bed  heater  to  deliver  hot  flow  to  a  stationary 
linear  cascade  at  a  selected  pressure  and  temperature.  In  the  experiments, 
test  times  were  0.8s  with  steady  flow  achieved  after  0.2s. 

The  injection  of  the  coolant  was  synchronized  with  the  hot  flow. 

Thus,  aerodynamic  and  heat  transfer  data  were  taken  under  steady  flow 


conditions  for  0.6s.  The  short  duration  of  the  test  means  that  the  blade 
surface  was  nearly  isothermal  during  the  test.  The  isothermal  condition 
approximates  actual  blade  operating  conditions  and  meets  the  design 
objective  of  avoiding  significant  temperature  gradients  within  the  blade. 
The  isothermal  wall  conditions  allow  a  logical  application  of  local  heat 

transfer  results  to  the  design  of  film  cooled  blades. 

( 42) 

Pressure,  temperature  and  blade  length  modeling  were  used  in  order 
to  scale  actual  gas  turbine  conditions  corresponding  to  a  turbine  inlet 
temperature  of  1800  K  (2750°F)  to  a  cascade  inlet  temperature  of  450  K 
(320°F).  The  temperatures  of  the  blade  surface  and  of  the  coolant  at 
injection  corresponded  to  two-thirds  and  one-half  of  the  inlet  absolute 
temperature,  respectively.  The  Reynolds  number  based  on  blade  chord  was 
5x10^.  The  turbulence  level  of  the  main  inlet  flow  was  10%.  The  details 
of  the  blowdown  facility  and  the  flow  modeling  are  given  in  Ref.  42. 

The  cascade  was  instrumented  with  total  pressure  and  temperature 
probes  located  at  a  distance  of  three  axial  blade  chords  upstream  of  the 
nine  blade  cascade  (Figure  48).  A  traversing  mechanism  carried  a  total 
pressure  directional  probe  at  mid-span  with  a  speed  of  2.2  cm  per  second. 
During  the  steady  flow  part  of  the  test  time,  the  probe  covered  a  distance 
larger  than  the  blade  wake  at  a  distance  from  the  trailing  edge  of  one- 
third  of  an  axial  chord. 

All  tests  were  performed  at  closely  identical  thermodynamic  conditions 
for  the  hot  flow  and  with  blade  surface  close  to  300  K  and  a  coolant 
temperature  of  225  K. 

The  blade  under  study  had  a  thick  trailing  edge  and  straight  suction 
side  (Figure  49)  and  was  previously  investigated  in  the  M.I.T.  blowdown 
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facility  by  Demuren  and  Hajjar,  et.  al.  for  both  off-  and  design 

conditions.  The  thick  blade  provided  space  for  the  cooling  plenums  while 

leaving  an  adequate  structure  for  a  thick  skin.  The  locations  of  the  line 

of  film  cooling  holes  (Figure  49)  were  selected  to  match  the  hot  spot 

regions  on  the  blade  surface.  Therefore,  five  single  rows  of  holes  were 

drilled  at  the  leading  edge.  At  the  transition  point  on  the  suction  side, 

at  the  mid-chord  on  the  pressure  side,  and  upstream  of  the  trailing  edge  on 

both  sides,  the  coolant  was  injected  through  a  double  line  of  holes.  The 

angles  of  injection  with  respect  to  the  blade  surface  were  of  20°,  based  on 

flat  plate  maximum  effectiveness  except  at  the  leading  edge  and  on  the 

suction  side  where,  because  of  manufacturing  constraints,  these  angles  were 

90°  and  43°,  respectively.  The  diameter  of  the  holes  was  1mm  and  the  hole 

interspace  was  equal  to  two  hole  diameters.  The  double  line  of  holes 

located  at  the  transition  point  corresponded  to  the  joint  of  minimum 

pressure.  The  double  line  of  holes  upstream  of  the  trailing  edge  was 

(45) 

located  using  the  results  of  Amana’s  study  on  trailing  edge  heat 
transfer . 

Among  the  nine  aluminum  blades  in  the  cascade,  the  three  at  the  center 
were  film  cooled  and  only  the  middle  one  was  instrumented  (Figure  3). 
Fourteen  surface  static  pressure  taps  and  eighteen  heat  transfer  gauges 
allowed  the  determination  of  aerodynamic  and  heat  transfer  characteristics. 

The  pressure  taps,  located  along  the  center  line  of  the  blade,  were 
connected  by  using  holes  drilled  spanwise  in  the  blade  skin,  to  pressure 
transducers . 

The  heat  flux  into  the  blade  surface  at  a  given  location  was  measured 
with  a  calorimetric  gauge  mounted  flush  with  the  surface.  The  heat 
transfer  gauge  consists  of  a  cylindrical  copper  slug,  0.24  cm  diameter  and 
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0.05  an  in  thickness,  insulated  from  the  blade  and  having  a  thermocouple 
attached  to  its  backside  to  measure  the  slug  temperature.  If  the  heat  flow 
into  the  blade  is  modeled  as  one-dimensional,  and  the  insulation  is  assumed 
to  be 

perfect,  then  a  heat  balance  for  the  gauge  is  given  by: 

q  =  h(T  -  T  )  =  d/dt  (mC  /A  T  )  (52) 

H  00  g  P  g 

or 

h  =  mC  /A  (T  -  T  )  dTg/dt  (53) 

P  ®  g 

By  measuring  the  local  slope  of  the  gauge  and  main  flow  stagnation 
temperature  under  steady  flow  conditions,  the  heat  transfer  coefficient  car. 
be  obtaintained  since  mC^/A  is  a  constant. 

The  heat  losses  from  the  gauge  are  also  minimized  by  designing  the 
insulation  thickness  to  have  a  longer  diffusion  time  than  the  test  time. 

The  Nusselt  numoer  is  calculated  as: 

Nu  =  hc/k  (54) 

where  k  is  based  on  the  mainstream  stagnation  temperature  and  the  blade 
chord,  c,  of  6.5  cm. 

The  study  of  the  wake  profiles  allowed  us  to  investigate  the  effects 
of  the  cooling  on  the  blade  profile  loss.  A  stagnation  pressure  loss 
coefficient,  5  ,  was  obtained  by  taking  the  ratio  of  the  stagnation 
pressure  drop  over  the  ideal  dynamic  pressure  present  in  the  flow  (in  the 
absence  of  any  loss  mechanism) ,  averaged  over  a  pitch  length  and  weighted 
by  the  mass  flow  as: 
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The  effect  of  cooling  on  the  heat  transfer  was  analyzed  with  the 
isothermal  effectiness  n  variations,  defined  in  Ref.  42  as: 
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where  q  =  heat  transfer  rate  to  the  blade  without  cooling 

o 

H 

qc  =  heat  transfer  rate  to  the  blade  with  cooling. 

For  each  of  the  seven  following  cooling  configurations,  six  blowing 
ratios  have  been  investigated  (m  =  (p  U > _/ (0  U >»  =  0;  0.3;  0.5;  0.7;  1;  1.2). 
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The  injection  eonfiguations  (Figure  49)  are: 

1)  leading  edge  injection,  plenum  A 

2)  suction  side  injection,  plenum  B 

3)  pressure  side  injection,  plenum  C 

4)  trailing  edge  injection,  plenum  D 

5)  leading  edge  and  suction  side  injection,  plenums  (A  +  b) 

6)  leading  edge  and  pressure  side  injection,  plenums  (A.  +  C) 

7)  full  injection,  plenums  (A  +  ,B  +  C  +  D) 

/ 

Nitrogen  refrigerated  at  the  desired  temperature  of  225  K  was  used  as 
a  coolant.  Nitrogen  was  cooled  by  passing  it  through  a  coil  immersed  in 
liquid  nitrogen  and  fed  into  one  of  the  four  tanks.  Each  plenum  (A,  B,  C, 
or  D)  had  its  independent  cooling  system  fed  from  one  of  the  four  tanks. 

The  flow  was  metered  into  each  individual  plenum.  The  coolant  was 
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introduced  from  top  and  bottom  in  each  plenum  in  order  to  maintain  low  flow 
velocities  and  attain  uniform  injection.  Each  plenum  was  also  instrumented 
with  its  own  thermocouple  and  pressure  gauge  so  that  the  flow  properties  of 
the  coolant  were  exactly  known  at  injection. 

5.3  Experimental  Results 

5.3.1  Pressure  Distribution 

The  static  pressure  distribution  around  the  blade  was  found  to  vary 

(47) 

slightly  with  the  coolant  mass  flow  injected  .  Figure  50  gives  the 
pressure  distribution  for  different  values  of  the  injection  parameter  (m) 
when  coolant  was  injected  through  all  the  plenums. 

The  blade  passage  cross  section  can  be  reduced  either  by  the  local  jet 
penetration  or  by  the  existence  of  a  thick  boundary  layer  fed  by  the 
coolant.  Consequently,  the  pressure  in  the  blade  passage  will  rise,  for  a 
constant  pressure  drop  across  the  cascade,  as  the  mass  flow  decreases  in 
the  film  cooled  blade  passages.  The  pressure  distribution  on  the  pressure 
side  near  the  leading  edge  was  also  found  to  vary  because  of  a  change  in 
incidence  associated  with  changes  of  the  effective  shape  of  the  leading 
edge. 

5.3-2  Effects  of  the  Leading  Edge  Injection 

On  the  pressure  side,  the  injection  of  the  coolant  causes  the 
transition  to  a  turbulent  boundary  layer  to  occur  earlier  than  for  the 
uncooled  blade  (around  s/L  =  0.2  instead  cl'  s/L  =  0.3).  The  roughness  of 
the  mid-chord  coolant  orifices  trips  (Figure  51)  the  boundary  layer  and  a 
high  heat  transfer  is  locally  recorded.  For  the  highest  coolant  ma^s  flow 
injected,  good  cooling  is  obtained  as  far  as  the  trailing  edge. 
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On  the  suction  side  (Figure  51)  at  low  blowing  value  (rn  =  0.3)  a  poor 
cooling  is  obtained  due  to  the  rapid  mixing  of  the  small  amount  of  injected 
coolant,  while  at  the  highest  blowing  rate  a  large  heat  transfer  reduction 
(up  to  50%)  is  measured  along  ♦•he  blade. 

The  wake  profiles  displav,  for  all  m,  a  shape  similar  to  the  uncooled 
blade  wake.  Changes  in  the  location  of  the  transition  point  due  to  coolant 
injection  have  a  negligible  effect  on  the  wake  and  losses  (Figures  62,  63). 


5.3.^  Effects  of  the  Suction  Side  Injection 

The  injection  of  coolant  on  the  suction  side  occurs  near  the  maximum 
suction  point  and  the  coolant  is  therefore  introduced  in  an  adverse 
pressure  gradient. 

The  measurements  of  the  wake  width  (Figure  52)  indicate  that  even  for 
the  small  values  of  the  blowing  factor  (m  =  0.3),  the  boundary  layer 
thickness  is  greatly  increased.  This  thick  boundary  layer  with  the  cooler 
gas  near  the  wall  leads  to  high  values  of  the  film  cooling  effectiveness  as 
far  as  the  trailing  edge  located  at  a  distance  equivalent  to  eighty  hole 
diameters  (Figure  53). 

The  depth  of  the  wakes  also  indicate  that  the  stagnation  pressure  loss 
is  substantially  increased  (Figures  62,  63). 

As  the  blowing  factor  increases,  the  wake  is  found  to  become  wider  and 
deeper  up  to  a  value  of  m  =  1.2,  where  the  wake  width  is  shown  to  be 
reduced  because  the  boundary  layer  energization  (by  the  cooling  jets) 
becomes  more  important. 
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5.3.5  Effects  of  the  Presure  Side  Injection 

On  the  pressure  side  (Figure  5*0,  downstream  of  the  injection  rt:ion, 
a  fairly  good  heat  transfer  reduction  is  obtained.  At  a  low  coolant  mass 
flow,  the  Nusselt  number  rises  sharply  behind  the  injection  location.  As 
the  coolant  mass  flow  increases,  better  cooling  is  achieved  but  the  rate  of 
heat  transfer  reduction  levels  off  beyond  m  =  0.7.  The  jet  penetration  is 
evidenced  at  m  =  1  and  1.2  by  the  minimum  of  heat  transfer  rate  observed  at 
the  second  gauge  indicating  the  zone  of  jets  reattachment.  Close  to  the 
injection  holes  an  optimum  blowing  factor  of  m  =  1.0  is  observed,  whereas 
beyond  the  second  gauge  the  thermal  effectiveness  increases  with  coolant 
flow  rate. 

5.3.6  Effects  of  the  Trailing  Edge  Injection 

The  instrumentation  of  this  section  of  the  blade  shows  a  quite 
remarkable  heat  transfer  reduction  (Figure  55).  The  isothermal 
effectiveness  reaches  98?. 

At  low  values  of  the  blowing  ratio  m,  the  wake  profile  deepens  and 
widens  a  little.  At  m  =  1.2,  the  wake  width  and  depth  are  reduced  because 
of  a  local  energization  of  the  boundary  layer.  Therefore,  the  total  loss 
coefficient  is  found  to  have  a  maximum  at  a  value  of  m  close  to  1.0. 

5.3.7  Effects  of  the  Combined  Injections 

The  effect  of  combined  leading  edge  and  suction  side  injection  around 
the  blade  nose  is  similar  to  the  sole  plenum  A  injection;  while  around  the 
plenum  B,  the  rate  at  which  the  heat  transfer  increases  downstream  of  the 
injection  is  observed  to  be  slightly  less  than  for  the  injection  through 
plenum  B  alone ,  (Figure  59).  The  wake  profiles  observed  in  this  case 


and  in  the  sole  plenum  B  injection  are  almost  identical.  However,  higher 
loss  coefficient  follows  (Figures  62,  63). 

The  combined  injections  of  the  leading  edge  and  the  pressure  side 
provide  a  heat  transfer  reduction  around  the  leading  edge  and  at  the  blade 
mid-chord  similar  to  those  observed  for  plenums  A  and  C  tested  individually. 
The  wake  profiles  follow  the  same  trend  and  the  higher  loss  coefficient 
varies  linearly  with  the  coolant  mass  flow  (Figures  62,  63). 

5.3.8  Effects  of  the  Full  Injection 

The  reduction  of  the  Nusselt  number  (Figure  56)  is  found  to  be  domi¬ 
nated  by  the  nearest  upstream  injection.  Around  the  trailing  edge  very 
good  cooling  is  observed  on  both  sides  and  overcooling  of  the  blade  occurs 
at  m  =  1.2.  The  distribution  of  the  reduction  of  heat  transfer  with  the 
cooling  indicates  a  fairly  judicious  selection  of  the  location  of  the 
coolant  injection.  The  wake  profiles  are  found  to  be  close  to  those 
observed  for  the  other  case  of  injection  on  the  suction  side.  The  total 
loss  coefficient  is  found  to  be  larger  than  in  any  other  configuration 
(Figures  62,  63). 

5.3.9  Variation  of  Isothermal  Effectiveness  with  Blowing  Ratio 
and  Comparison  with  Flat  Plate  Correlations 

The  variation  of  the  isothermal  effectiveness  n  with  the  blowing 
parameter  m  is  shown  in  Figures  57  to  60.  Around  the  leading  edge 
(plenum  A)  on  both  sides  (Figure  57),  n  varies  linearly  with  m.  However,  a 
higher  effectiveness  is  achieved  on  the  pressure  side  because  of  a  larger 
number  of  coolant  orifices. 
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On  the  suction  side  (plenum  B)  (Figure  58),  the  effectiveness  reaches 
a  maximum  value  which,  close  to  the  injection  (s/D  <  12),  is  obtained  for 
m  =  0.7  and  for  m  =  1.,  further  downstream. 

On  the  pressure  side  (plenum  C)  (Figure  58),  the  rate  of  change  of  t.he 
effectiveness  with  m  varies  significantly  from  location  to  location. 

Around  the  coolant  outlets,  n  reaches  a  maximum  for  m  beyond  0.7.  Further 
downstream  (s/D  >  10)  n  increases  at  a  smaller  rate  and  no  maximum  value  is 
observed. 


For  combined  and  full  injection  (Figures  59,  60),  the  curves  display  a 
similar  trend  in  the  variation  of  n  with  m.  As  above,  an  upstream  injec¬ 
tion  does  not  significantly  influence  the  effects  of  downstream  injection. 

An  attempt  was  made  to  compare  the  recent  data  with  the  correlation 
introduced  in  Ref.  (42,43),  for*  isothermal  effectiveness  measured  dowstream 
of  the  line  of  holes  and  is  shown  in  Figure  61.  It  displays  the  variation 
of  the  isothermal  effectiveness  versus  A^  as  defined  (Figure  61) 

as  A  =  k  “1#35  8~1,35  Re  ~0’25  (Cos2a)-1,35  s/D  (57) 

u  u  o 

The  law  for  injection  v.-  •  jgh  a  double  row  of  holes  on  flat  plate  without 
pressure  gradient  has  .  correlated  as  follows: 

n  =  1.46  AQ  eXp(',25)  (58) 

This  is  represented  by  a  straight  line  on  Figure  61.  The  data  obtained  on 
the  suction  side  lies  slightly  above  it.  As  commonly  found  by  experimental 
researches,  injection  on  a  convex  wall  improves  the  cooling. 

On  the  blade  pressure  side  (plenum  C)  the  lower  effectiveness  confirms 
that  the  coolant  undergoes  a  larger  mixing  than  when  injected  through  a 
slot  on  a  flat  plate  in  the  absence  of  a  pressure  gradient. 

Similar  discrepancies  occur  when  the  results  of  a  different  corre¬ 


lation  which  uses  the  concept  of  an  equivalent  width  for  injection  through 
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holes v  are  compared  with  the  data.  These  discrepancies  are  not  un¬ 
expected  since  individual  jets  issued  from  the  holes  penetrate  deeper  into 
the  hot  flow  on  the  suction  than  on  the  pressure  side  for  the  same  value  of 
the  blowing  ratio  Cm)  because  the  variation  of  pressure  between  the  two 

« 

sides  leads  to  a  difference  in  momentum  normal  to  the  flow  for  the  same 
value  of  m. 

( -so ) 

A  calculation  of  the  individual  jet  penetration  indicates  that  the 
jet  center  penetrates  up  to  1.2  hole  diameters  on  the  suction  side  as  com¬ 
pared  to  0.76  hole  diameters  on  the  pressure  side  for  the  same  value  of  the 
blowing  ratio  m  =  1.0. 

The  depth  of  penetration  of  the  jet  and  the  local  dynamic  pressure 
influence  the  drag  of  the  jets  on  the  hot  flow,  which  will  therefore  in¬ 
crease  with  the  depth  of  penetration  or  with  the  blowing  ratio. 

This  mechanism  appears  to  describe  the  quasi-linear  dependence  of  the 
profile  losses  on  the  blowing  pressure  observed  when  films  are  injected 
either  from  the  pressure  or  suction  side.  For  injection  upstream  of  the 
trailing  edge,  a  similar  mechanism  is  likely  to  dominate  for  small  values 
of  the  blowing  ratio.  For  larger  values  of  m  the  film  energizes  the  wake 
and  reduces  the  losses. 

The  depth  o’  penetration  of  the  individual  jets  influences  the 
thickness  of  the  mixing  layer  formed  when  the  jets  are  in  contact  with  each 
other  and  with  the  wall.  The  thickness  of  this  layer  is  larger  on  the 
suction  side  than  on  the  pressure  side  where  the  layer  is  quickly 
accelerated  and  subjected  to  rapid  mixing  under  the  influence  of  a  strong 
negative  pressure  gradient. 
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5.4  Conclusions 

1 )  Static  pressure  distribution  around  the  blades  varies  only 
slightly  with  the  coolant  mass  flow  injected. 

.  2)  The  leading  edge  injection  provides  at  high  injection  (m  =  1.2) 

effective  cooling  along  the  blade  wall  with  low  stagnation  pressure  loss. 
The  dependence  of  the  profile  losses  is  found  to  be  non-linear  with  the 
mass  flow  rate  of  coolant,  because  they  are  influenced  by  mixing,  motion  of 
the  transition  point  and  changes  in  the  location  of  the  stagnation  point. 

3)  Injection  on  the  suction  side  leads  to  important  heat  transfer 
reduction,  but  also  to  the  largest  loss  among  the  cases  tested.  The 
profile  losses  are  found  to  be  linear  to  the  rate  of  coolant  injected. 

4)  Injection  on  the  pressure  side  provides  good  cooling  close  to  the 
injection  zone.  Further  downstream  strong  mixing  dominates,  while  losses 
are  linearly  dependent  upon  the  cooling  mass  flow  rate. 

5)  Trailing  edge  cooling  is  achieved  effectively  on  both  sides.  The 
profile  losses  are  found  to  vary  linearly  with  the  coolant  mass  flow  rate 
up  to  a  value  of  m  =  1.0  when  the  profile  losses  are  reduced  due  to  the 
energization  of  the  wake  by  the  coolant. 

6)  When  combined  injection  is  used,  the  heat  transfer  reduction 
achieved  downstream  of  each  coolant  outlet  is  mostly  determined  by  the 
effect  of  the  closest  upstream  injection.  Full  injection  leads  to 
important  losses  in  stagnation  pressure;  e.g.,  for  a  value  of  the  blowii.g 
parameter  m  a  1.0,  the  overall  profile  losses  are  eighty  percent  larger 
than  for  the  case  with  no  cooling. 
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5.5  Nomenclature 


A  =  heat  transfer  slug  area  open  to  the  main  stream 
c  =  blade  chord 

Cp  s  specific  heat 

D  =  hole  diameter 

h  =  convective  heat  transfer  coefficient 

k  =  thermal  conductivity 

2  2  1/2 

kQ  =  momentum  ratio  parameter  =  (PU  )c/(pU  )ra 

L  =  total  length  along  the  blade  from  L.E  to  trailing  edge 

m  =  blowing  ratio  m  =  (PU)c/(oU)0# 

m  =  coolant  mass  flow 

c 

m  =  main  stream  mass  flow 

00 

Nu  =  Nusselt  number 

P  s  static  pressure 

qQ  =  heat  flux  to  the  blade  without  cooling 

qc  =  heat  flux  to  the  bl’ade  with  cooling 

Rec  =  Reynolds  number  based  on  coolant  conditions  and  injection  slot  or 
hole  dimension 

s  s  curvilinear  coordinate  along  the  blade  from  the  leading  edge 
T  =  main  stream  total  temperature 

Tc  =  coolant  temperature 

T  =  heat  transfer  slug  temperature 

ft 

U  =  velocity 

Subscript 


0  =  stagnation  condition 

1  =  position  upstream  of  the  cascade 

2  =  position  downstream  of  the  cascade 

Greek 

a  =  angle  of  injection 

=  stagnation  pressure  loss  coefficient 
p  =  density 

3  =VT- 

Other 

Aq  s  empirical  correlation  coefficient  for  isothermal  effectiveness 
on  o  flat  plate 
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Fig,  2q,  Jet  Cross  Section  and  Curvilinear  Coordinates 


Fig.  25,  Sections  of  the  Jet 
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Fig.  27,  Curve  ^ Along  Which  Kelvin's  Circulation 
Theorem  is  Applied 
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FlG  28,  CONTRAROTATIVE  VORTICES  IN  ANY 
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Fig. 30,  Reattachment  of  the  Jets 
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Fig.  32,  Velocity  Induced  on  a  Vortex  in  the  Cross  Section 
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Fig.  35.  Concentration  Profile 


Fig,  36,  Evolution  of  the  Velocity  Profile 


Fig.  39.  Trajectory  of  a  Vortex  in  the  Layer 
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Fig. 40.  Concentration  Profile  in  the  Layer 
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